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Abstract 

In this paper we study an initial boundary value problem for Za- 
kharov's equations, describing the space propagation of a laser beam 
entering in a plasma. We prove a strong instability result and prove 
that the mathematical problem is ill-posed in Sobolev spaces. We also 
show that it is well posed in spaces of analytic functions. Several con- 
sequences for the physical consistency of the model are discussed. 
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1 Introduction 

Zakharov's equations [12] model electronic plasma waves, describing the 
coupling between the slowly varying envelope of the electric field E and the 
low-frequency variation of the density of the ions n. A commonly used form 
of the equations reads 



(1.1) 



i(d t + ^d z )A + -^—A X A = ^5nA, 
ujq Zuiq 2noWo 



where loq is the frequency of the laser, ko its wave number and uj pe the 
plasma electronic frequency; they are linked by the dispersion relation Uq = 
u!p e + k^c 2 where c is the speed of light; uq the mean density of the plasma, 
rrii is the mass of the ions and c s the sound velocity in the plasma. The space 
variables are (z,x), z £ M and x E M. 2 ; z is the direction of propagation of 
the laser beam and x are the directions transversal to the propagation. In 
this model, the transversal dispersion is neglected. 

Reduced to a dimensionless form (see in Section 6) , the equations become 



(1.2) 



i{ed t + d z )E + A X E = nE, 
(8 2 - A x )n = A X \E\ 2 



See [9] or [TO] for the introduction of this kind of models for numerical 
simulation. 
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The local in time Cauchy problem for (jl.2p is now well understood: in 
[6], it is proved that it is well posed, locally in time, for data in suitable 
Sobolev spaces. This extends previous results for the classical Zakharov 
system, where transversal dispersion is taken into account, that is when A x 
is replaced by A( z ,x) ( see PHI El E] and references therein). However, the 
system f)1.2|) is quite different from the classical Zakharov system since the 
Cauchy problem for periodic data exhibits strong instabilities of Hadamard's 
type ([!]). Note that periodic data are considered in the quoted paper as a 
model for data which do not vanish at infinity. We refer to [I] for a discussion 
of the physical relevance of the different frameworks. In particular, we recall 
in Section 6, why the periodic context is well adpated to physical situations, 
where the envelop of the beam has rapid oscillations (speckles). 

In this paper we consider a boundary value problem for (jl.2j) which mod- 
els the the propagation of a laser beam entering the plasma at an interface 
{z = 0}. This approach is very common in physics where people is actually 
more interested in describing the propagation in space rather than in time, 
i.e. considering z as the propagation variable. Indeed, this is an underlying 
idea in many of the paraxial approximations, like the one which yields the 
Schrodinger equation in (jl.2p . This idea is also very common in numerical 
simulations. In this case, the system (j 1 . 2 h is considered in the half space 
{z > 0}, for positive times {t > 0} together with initial-boundary conditions 

(1.3) n\ t=0 = 0, d t n\ t=0 = 0, E z=0 = E , E t=0 = 0. 

This approach is also very natural when the parameter e is small. The 
solutions are expected to vanish for t < z, by finite speed of propagation. 
Changing t to t — ez, the system reads; 



id z E + A X E = nE, 
{d 2 t -A x )n = A X \E\ 



(1.4) 



(1.5) 



We look for solutions U = (E, n) of (jl.4p (jl.5p which are periodic in x 
with period L. We denote by T the corresponding torus (W/2-kL) 2 . As 
recalled in Section 6, periodicity is somewhat natural when the envelop of 
the laser beam has a transversal structure which involves length-scales that 
are large with respect to the length-scale of the laser but small compared 
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to the width of the beam (see [3]). Periodicity is also natural for the use of 
spectral methods in numerical simulations. 

The main result of this paper is a strong instability result for (jl.4p . in 
the spirit of @]. It says that the boundary value problem (|1.4p (jl.5p is ill- 
posed in Sobolev spaces. Prior to this result, we make a detailed analysis of 
amplification properties of the linearized equations : at a space frequency of 
length k, the amplification is exponential in e cv ^ . This implies a Hadamard's 
type instability and the ill-posedness in Sobolev spaces, as well a the local 
well-posedness in spaces of real analytic functions. For physical interpre- 
tations the frequencies must be considered as bounded while for numerical 
applications they are filtered by one way or another. Thus is is important 
to give a numerical value to the amplification rate and to be evaluate the 
amplification effect on definite values of the fields. This is done in Section 6. 

Let us describe now the main contents of the paper. We first consider 
the linearized equations. Consider here a constant solution 



which satisfies (jl.4j) (II. 5ft with boundary data Eq = E. The homogeneous 
linearized equations near U_, at the frequency £ £ M 2 read: 

( id z u — k 2 u — En = 0, 



where k := \£\ and u [resp. v] [resp n] denote the the Fourier coefficient at 
the frequency £ of the variation of the field E [resp. E] [resp of the density]. 
They are supplemented with initial boundary boundary conditions: 



The equations (11. 7p (|1.8p form a well posed hyperbolic Goursat problem (see 
PQ). Our concern is to understand its behavior for large k. The fundamental 
solution is studied in details in Section 2, where we also construct blowing 
up solutions: 

Theorem 1.1. There are initial data \uq\ < 1 and vq = 0, such that for 
k large, t > 4 and p := yj2k\E\ 2 zt large, the solution of the homogeneous 
equation (jl.7p satisfies 



(1.6) 



(1.7) 




(1.8) 



n\ t=0 = 0, d t n\ t=0 = 0, u z=0 = u , v z=0 = v . 



(1.9) 



n(t,z)\ > kp 2e p . 
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The exponential amplification in e c * is the signal of a strong instability. 
We construct solutions on domains 

(1.10) O = {(t, z, x) G [0, T] x [0, Z] x T 2 ; zt < 5} 
and prove that they blow up on the part of the boundary 

(1.11) r = {(t,z,x) £ n ; zt = 5} 
which is not empty if TZ > 5. 

Theorem 1.2. For all s, T > L > andE_ ^ 0, there are sequences 5 k and 
Z k and families of solutions Up. = U_~\~ i e ki n k) °f (jl-4|) . in C ^^ H S {T)) 
such that 



(1-12) l|6A;|2=ol| J H' s ([0,T]xT 2 ) ~ ► °> 

(1.13) ^fc^O, 5 fc ^0, TZ k >5 k , 

(1.14) sup ||n fc (t,2:)|| L 2( T ) oo. 

(t,z)er k 



This theorem is proved in Section 3, with technical details postponed 
to Section 4. This nonlinear instability result is pretty strong: not only 
the amplification ||e||i2/||e| z= o||H s is arbitrarily large, in arbitrarily small 
distance Z, with arbitrary loss of derivatives s, but there is an effective blow 
up of the norm of n k . 

This analysis reveals the importance of the amplification factor, p := 
y / 2k\E\ 2 zt, and indeed there is a good uniform stability for a filtered system 
at frequencies |£| < k, on the domain {(t, z, x) G [0, T] x [0, Z] x T 2 } provided 
that 

(1.15) kZT\E\ 2 < 1. 

Instead of filtering the frequencies, another mathematical approach is 
to counterbalance the amplification by an exponential decay of the Fourier 
coefficients. This means that one works in spaces of real analytic functions. 
In this framework, we prove in Section 5 a local existence theorem by an easy 
adaptation of proof of the existence of analytic solutions to the hyperbolic 
Goursat problem (see |llj). The interesting point is that the length Z of 
propagation satisfies an estimate which is very similar to (|1.15p : 

(1.16) RZTWE^f^ « 1. 

where E^j is a space of analytic functions for the boundary data on [0, T] x 
T 2 and iT 1 measures the width of the complex domain where Eq can be 
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extended. Indeed, for Eo(t,x) = eo(t)e l ^' x with |£| = k, and R k, k, there 
holds H-EoIIeh <E_ = ||e||ioo and (|1.16p is equivalent to (|1.16p . 

This shows that (|1.15p or (|1.16p can be seen as stability criteria for the 
Goursat problem for Zakharov equations (jl.2p . 

Section 6 is devoted to a qualitative discussion of the results. We discuss 
several points. 

- Taking £ to be a constant is not physically realistic : it would mean 
that the envelop of the mean electric field has a jump. The case where E 
is a smooth function of time will be briefly discussed. It yields additional 
technical difficulties but does not change qualitatively the results. 

- The boundary data in Theorem 11.11 are very particular. Thus, it 
is important to understand better how general solutions behave. Actually, 
the rate of amplification also depends on the time frequencies. The under- 
lying phenomenon is a resonance between plasma waves governed by the 
Schrodinger equation and acoustic waves for n. In other words, only the 
acoustic oscillations e - lujt+t € x with lo 2 = |£| 2 which present in the boundary 
data Eq are effectively amplified at the given exponential rate. Mathemat- 
ically, in general, these acoustic frequencies have a nonvanishing amplitude 
because the signals exactly vanish in the past and thus their Laplace-Fourier 
transform has no lacuna. Physically, the acoustic frequencies, even when 
absent from the main scene, can be present in a background noise. By a 
standard plane wave analysis, we will also give an amplification rate for os- 
cillations which are not exactly acoustic, but this is not a correct approach 
for the Gourset problem. 

- The blow up in Theorem 11.21 is totally unphysical, since n is a vari- 
ation of density and thus must remain bounded. However, for physical 
interpretations, one must keep in mind that that (11. ip is only a model which 
has a limited range of validity. In particular, it is tacitly assumed that the 
variations of the ion density are not too large and that the paraxial approx- 
imation for the envelop is valid. Moreover, the physical frequencies k are 
bounded. One crucial question is to know wether the factor p = ^j2k\E\ 2 zt 
is small or large. In Section 6, we will give standard physical data for laser- 
plasma propagation showing that this factor can be large (10 to 10 2 ) for 
not very intense fields E of order 10 9 or 10 10 Wm" 1 . In this case, that the 
amplification e p ranges from 10 4 to 10 8 . Thus, the acoustic boundary oscil- 
lations can be ignored only if their relative value is <C e _p , which is much 
beyond the usual admissible errors. This seems to indicate that the model, 
as is it, is not well adapted to the propagation of intense laser beams. 

- For numerical simulations, the analysis shows that increasing the 
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number of Fourier modes, which is natural to improve the accuracy of com- 
putations, may introduce strong instabilities. Of course, one can eliminate 
most of them by filtering out the bad acoustic oscillations, but then the 
question is the relevance of the computations with respect to the model. 



2 The linear instability 

Consider the linearized equation from (jl.4p around (E_, 0) : 
I id z e + A x e — En = /, 



(2.1) 



{d 2 t 



A x )n - A x (Ee + Be) = A x h 



(2.2) 



n\ t=0 = 0, d t n\ t=0 = 0, e z=0 = e . 



Multiplying e by a constant phase factor e ld , there is no restriction in as- 
suming that E is real. Taking u = e and v = e as (independent) unknowns, 
the system reads 



(2.3) 



' id z u + A x u — En = /, 
id z v — A x v + En = g, 
{ (d t 2 - A x )n - EA x {u + v) 



A x h, 



with g = —f. 

Performing a Fourier series expansion in x (or a Fourier transform), and 
still denoting by u, v and n the Fourier coefficients, the equations at the 
frequency (6K 2 read with k := |£|: 



(2.4) 



id z u 

id z v 
>2 , ,:2 



k 2 u 



En = f, 



k 2 v + En = g, 



-k 2 h, 



{ (d'f + k z )n + k 2 E(u + v) 
together with initial boundary boundary conditions: 
(2.5) ri| i=0 = 0, d t n\ t=0 = 0, u z=0 = u , v z=0 



For L-periodic functions, the frequencies ^ G ^Z 2 . 

The Goursat problem (|2.4p (|2.5p is well posed (jlj). The main purpose 
of this section is to prove estimates for the fundamental solutions (Proposi- 
tions [23] and [2JJ] below), and give an example of a solution of the homoge- 
neous equation which is amplified at indicated rate indicated in Theorem ll.il 
(see also Theorem 12.81 below for a more precise statement). 
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2.1 The fundamental solution 

Extend the functions by for t < and perform a Fourier-Laplace transform 
in time; this amounts to replace dt by iC, with £ lying in the lower half 
plane {Im£ < 0}: the third equation in (|2.4|) and the homogeneous initial 
conditions for n imply that 

k 2 

(2.6) h = k2 (E(u + v)+h). 

We denote here by (f>(z,C) the Fourier-Laplace transform of (j)(z,t). We end 
up with the system 

(2.7) d z U = iAU + F, U {z=0 = U 
for 

with 

(2.9) A(C) 1 - k2 ~ a k2lKl 



a a + k 2 J ' C 2 - k 2 ' 

Note that j4(C) is bounded and holomorphic for {Im£ < —7} for all 
7 > and therefore, by inverse Laplace transform: 

Lemma 2.1. The solution of the homogeneous system (|2.4p with initial- 
boundary values (12. 5p where Uq = t (uQ,vo) G Co°(IR + ) on {z = 0}, is 

U(t,z) = ±- f e^^MOdC, 

n(t,z) = ^[ l,'M%„(()Ad(, 
2vr 7 M _ i7 C 2 - k 2 

where I := (1,1) and 7 is any positive real number. 

Note that the integrals above are convergent, as a consequence of the 
estimates given below. 

The fundamental solution is therefore linked to the distribution 
(2.10) £(t,z) = — I e %+^«)) d(, 

27T 7 K _j 7 
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where the integral is taken in the sense of an inverse Laplace transform. 
More generally, we are led to consider integrals of the form 



(2.11) I(t,z) = ±- f e^^piOdC, 

for rational functions p. 

Note that a and A and thus e*(*^ +2tA ^) are holomorphic in C\{— k, +k}. 
In addition, a(£) = O ( | C I 2 ) a t infinity, and thus 

(2.12) A(Q -A(oo) = 0{\C\~ 2 ), e izA{() _ e izA{oc) =0 (\C\- 2 ) 
where 

(2.13) A(oo) = 



-k 2 
k 2 



This implies the following 



Lemma 2.2. If p is a rational function with poles in {Im£ > 0} ; and 
p(C) = C C I CI 2 ) a t infinity, then the integral in (12. 8j) is absolutely convergent 
for all t G 1R. It vanishes for t < and for t > 0, 

(2.14) X (M) = ^ yV C+izA(C) KC) d( 

where T is any simple contour oriented positively winding around the poles 
of p and a. 

Proof. The integral is clearly convergent and letting 7 tend to +00 implies 
that it vanishes when t < 0. When t > 0, the integral over a large half circle 
in the upper half space tends to as the radius tends to infinity, allowing 
to close the integration path. □ 

Similarly, we can split £ into two parts: 
S(t, z) = — ! e^ + ^ (oo) ) d( + — [ e ^(e izA ^ - e izA ^) d(. 

Thanks to (|2.12|) . the second integral can be deformed to an integral on a 
closed contour T, on which the integral of the entire function e lt ( e lzA {° ) 
vanishes. Therefore, 
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Lemma 2.3. The distribution £ defined in ()2.10j) is equal to 

(2.15) £(t,z)=e izA ^ S t=0 +£ (t,z) 
where £$ vanishes for t < and is equal to 

(2.16) £ (t,z) = — [ e lK+lzAk ^ d( 

2vr J r 

when t > 0, where T is any simple contour oriented positively surrounding 
the poles of a. 



2.2 Estimates 



We give sharp upper bounds for the contour integrals in (|2.14p and (|2.16p . 
The matrix A is traceless, therefore A 2 = — (det A)Id = k 2 (k 2 + 2a)Id and 



JzA 



cos(kz\)ld + 



sin(kzX) 
kX 



A 



where 
(2.17) 



A = \/k 2 + 2a. 



The choice of the square root is irrelevant in the expression above of e tzA . 
Therefore, 



(2.18) 



.1 it^ikXz ( 2 +k + I 
1 fc _ A 
A fc 



A _ k 

k i A t 

r> A fe 

fc A 



_i_ ^ iti^+ikXz 

+ 4 



9 A k 

k X 
X _ k 
k X 



k _ X 
X k 

2 + I + ! 



This reveals the phases 

(2.19) <P = t£- zk\, 



if) = t( + zkX 



Lemma 2.4. For k > max{l, 4|i?| 2 } ; choosing the principal determination 
of the square root in (|2.17p defines X(k + w) as an holomorphic function of 
w for 1 < < ifc. Moreover, for all t > 0, z > and w in this annulus, 
there holds 



ip(k + w) — (kt - k 2 z + tw 



kz\E\ 
2w 
kz\E\ 2 



ip(k +w)- (kt + k 2 z + tw+ ^— 1 ) < ^(||^;| 2 + l\E\ 4 ). 
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Proof. For |to| < ^k, 

a{k + w) 

Thus 



k 2 \E\ 2 



k\E\ 2 \E\ 
w(w + 2k) ~ 2w 4(1 + 



2k- 



k\E\- 



2w 



< - \E\ 



If in addition \w\ > 1 and k > max{l, 4|£| 2 }, we have 



< \E\ 2 , 



2a 

k 2 



2\E\ 2 1 
< — 1 < -. 
~ k ~ 2 



Thus we can choose A to be the principal determination of the square root 
of k 2 + 2a and 



(2.20) 



a 2 „,2a. 



k\ = k 2 + a + -2) 



A; 2 K k 2 ' 

where G is holomorphic on the unit disc. Substituting, we find that 



(2.21) 



\k\-k 



2 k\E\ 2 > - 1 



2w 



<o\E\' + \E\" sup |G(C)|. 
6 Kl<l 



The supremum is \G(— k)\ < i and the lemma follows. 



□ 



Proposition 2.5. T/ien there is a constant C, such that for all E, k, t and 

z satisfying 



(2.22) 



1. 



\E\ < 1, k>m&x{l,4\E\ z }, ie[0,l], z G [0, -], 



i/ie contour integral in (|2.16p satisfies : 



\£o(t,z)\ < C 



r & 



with p = y/2k\E\ 2 zt and r = 1 + k ^ p z ■ 

Proof, a) We choose T to be the union of the circles 

ri = {|C-fe|=r}, T 2 = {\C + k\=r}, 
with r 6 [1, T^k] to be chosen depending on t and z. 
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By flgjZH and P^2j> - we have the following bounds: 



Furthermore 



ft; 

A" 1 



< 



2/cr 



+ 



A: 2 



< 



El 

ft; - 4' 



|£| 2 |£fz 2|e;| 2 i 

- 2|A|r + ftr|A| ~ k ~ 2" 



Thus 



A fc 
I + A 



<3, 



A ft; 
ft? ~ A 



< 



3|£| 



k 



A fc 
fc ~ A 



< 



3\E\ 



k 



Therefore, all the entries of the two matrices present in (|2,18p are bounded 
by 5. 

b) If 

(2.23) k\E\ 2 zt > 2, 

then we choose 



(2.24) 



kz\E\ 2 
2t 



By (I2T231) and fl!T22"1) . r G [1, \k}. For C = k + re* e G r*i, Lemma E3 implies 
that the imaginary parts of the phases satisfy 



(2.25) 
(2.26) 



limy) (A; + re ) - psin6>| < z\E\ , 



\lmip(k + re ie \ < z\E 



Therefore, the integral over T± contributes to £q to a matrix whose entries 
are bounded by 



1 

2^ 



2~ 



1 + e 



-psm 







? ) r dO. 



By symmetry, the integral over T2 is estimated similarly, and therefore, for 
p > 2, the entries of £q are bounded by 

»2tt 



5r 

47T 



, Z |E|2 



The integral is bounded by 

-it/2 



37T + 



e pcosy ^ = 37r + 2e^ 



-7T/2 



< 3tt + 2e p 



1 + e" 



tt/4 

e" 

-7I-/4 

\/2/2 

i 

-V2/2 



-psm 



? ) d0. 



-2psin y 



-2pt 2 



dy 



dt 



< 37T + 
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Therefore, for p > 2, the entries of £q are bounded by 



i - i - /.-- i\E\ 2 kz /3 e p \ \E\ 2 kz e p 

(2.27) \S {t,z)\<he z ^\ «/J=L_ (- + __) <5ViA 



t \4 l^/ixp) V t y/np 



c) If 



(2.28) -p 2 = k\E\ 2 zt < 2, 

2 

we choose 

(2.29) r = max{l,k\E\ 2 z} £[l,^k}. 

Consider first the integral on T±. The entries of the corresponding matrix 
The imaginary part of the phases are bounded by 

\E\ 2 

llm J < tr + zk ] -=^- < max{t, tk\E\ 2 z} + 1 < 2 + p. 
Therefore, the integral are bounded by 

^-e 2+ o < 10 ef. 
4 

The analysis of the integral over T2 is similar, and combining the estimates 
above, one obtains (|2.19p □ 

Proposition 2.6. Denote by E\ [resp. £2 ] the contour integral in (12. lip 
with p = £2_ k 2 [resp. p = ^_ k -iyi ]■ There is a constant C, such that for 
\E\ < I, k > max{l,4|£;| 2 } ; t € [0, 1] and 2\E\ 2 z < 1 there holds 

(2.30) \£i{t,z)\ < C- k "" 



(2.31) MM) I <C lr "' 



Proof. For £1 we proceed as for £0, noticing that |^rp| ~ f on the inte- 
gration path. 
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For (|2.3ip . we make the same choice of radius r when k\E\ 2 zt > 2. 
fc 4 I 



, 4 , 2 

Noticing that | lr -±_ h -M | ~ on the integration path yields the estimate 



|£ 2 (M)| <Ci 



t k 2 eP 



\E\ 2 kz y/T+~p' 



When k\E\ 2 zt < 2, the length T is 0(r), we maximize r in the range 
[1, 7}k] instead of minimizing it as we did in (|2.29p . so that the imaginary 
parts of the phases remain bounded. Namely, we choose 

f = minj — , — k\ 
x t ' 2 s 

and get 

k 2 

\£ 2 (t,z)\ < C-. 

in this zone. Combining the estimates yields (|2.31j) □ 

2.3 Integral formulas for the solutions of (12.41) (12.51) 

At the given frequency k, the system (|2.4p (|2.5|) is a classical hyperbolic 
Goursat problem in dimension 2, and therefore has a unique smooth solution 
for smooth data (see pQ), which respects the causality principle in time. We 
will use the following description: 

Proposition 2.7. For F = \f,g) and h in C°°([0,T] x [0, Z]) and the 
problem (|2.4p (|2.5p with Uq = *(uo,«o) = 0, has a unique solution U = 
\u,v), n in C°°([0,T] x [0,Z]) and 



U(t,z)=i e^ z - z '^ A ^F(t,z')dz' 
Jo 

(2.32) + / / £ {t!,z')F(t-t!,z- z')dz'dt' 

JoJo 

+ ^2 £l(*'> z ') (7) ^* ~ *'' z ~ z ') dz ' dt '' 



(2.33) 



n(t,z)= / / (— i, i) £i(t' , z')F(t — t' , z — z')dz' dt 1 
J0J0 

+ i £2(t' , z')h(t — t' , z — z')dz' dt' . 
J0J0 
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2.4 Exponentially growing solutions 

We show that the rate of amplification e p observed in Propositions 12.51 and 
is sharp. Consider the solution of 



' id z u + A x u — En = 0, 
(2.34) < id z v - A x v + En = 0, 

, (df- A x )n-EA x (u + v) = 0. 

with initial-boundary conditions 

o n s'm(kt) i£ 
n\ t=0 = d t n\ t=0 = 0, u\ z=Q = — ; e ? 



(2.35) 

This amounts to solve 
(2.36) < 



k 



' id z u — k 2 u — En = 0, 
d z v + k 2 v + En = 0, 



0. 



k (d t 2 - A x )n- Ek 2 (u + v 
with k = |£|, together with initial-boundary conditions 

(2.37) 



v\z=o = 0. 



ri| 4=0 = 9t^|t=o = 0, n^ =0 = — - — , v\ z=0 = 0. 



Because the Fourier-Laplace transform of l{t>o} — j: 
is 



— ' is fcTT^-) the solution 



(2.38) 

(2.39) 
where 



U(t,z) 
n(t, z) 



u(t,z)\ = J_ 

v(t,z)J 2i7T 



_i(tC+*A fc (0) 
-17 K S 



J), ^=(i,i)- 



Theorem 2.8. There are constants C > c > suc/t i/iat /or k > 2, 

te [o,i], ze [o,i] ; 

1 e p 
£/(i,z) <C--= 

KV1 + /0 

|n(* ) i)|<C(~ + -^-)-^=. 
1 1 x k I + p' yjl + p 
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Moreover, for p > 1, t > \ and k large enough, 
(2.40) \n k (t,z)\>cp-lef(l + 0(p- v 



Proof. One can deform the integration path R — to a contour T which is 
the union of the circles 

r 1 = {|c-fc|=r}, r 2 = {|c + fc|=r} 

with r G [1, The upper bounds follow from Proposition 12.61 

We now concentrate on the proof of the lower bound for n assuming that 
p = y / 2\E\ 2 zt > 2. We choose the radius r = £ as in (pT24"jh By (|27[8jh for 
j = 1,2, the integral over Tj contributes to 

Nj(t, z) =— I e*-***— ^ 



2 m } v . 2A (( 2 -k 2 ) 2 C 

If itc-izkx^ + k k 2 E d( 
2mj T . 2A (( 2 -k 2 ) 2 C " 

Consider first N\. We evaluate the phases (p = t( — zkX and if) = tQ + zk\ 
on T 1 as in ([235]) ([2T26]) . For C = A; + re ie we have 

(^(A + re* 61 ) = kt - k 2 z + 2ip sin 9 - a k (r, 9), 
ip{k + re ie ) = kt + k 2 z + 2p cos 6> + a k (r, 9). 

where the a k are smooth functions of 9 G T, uniformly bounded as well as 
their derivatives when 1 < r < ^k. 

Note that the imaginary part of tjj is bounded, while by (|2.2ip 



(2.41) 



X-k k 2 \E\ 2 



2A (C 2 - k 2 f 



, C|^| 2 C|£| : 



on Ti. Thus the corresponding term contributes to 0(k 1 ) terms in N\. It 
remains to evaluate 

We have the following expansion 

X + k k 2 E e~ 2ie 



2A (C 2 - k 2 ) 2 r 2 



Pk{r,9) 
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where the pk are smooth functions of 9 G T, uniformly bounded as well as 
their derivatives when 1 < r < and such that, uniformly on this domain 

E 

lim Pk(r,0) = — 

k— »oo Z 

Substituting, we see that 

i(kt—k 2 z) r2ir 

N x (z,t) = =— / e - 2psine Pk (r,e)e- lzr7k ^- 2t0 d9. 

The stationary phase theorem to implies that 

Nx(z,t) = e i ( fct -* a «) -^ = ( Pfc ( r ,_ 7r /2) e -^('-.-' r /2)+i 7 r + 0(p -i)) 

Therefore, for large enough, 

|JVi(f,*)| > J|^^(l + 0(p-i)) +0(1). 

b) The analysis of the contribution of the integral over T2 is similar, 
except that the phase with large imaginary part is now ip. Noticing that the 
corresponding amplitude is bounded by (|2.41|) . we see that 



r 2 y/p 



\N 2 (t,z)\ < -— + 0(1) 



and because i w & < ^, the estimate (|2.40p follows for fc and p large 
enough. □ 



3 The nonlinear instability 

3.1 The method of proof 

We fix a constant field E 7^ 0, noticing that (E, 0) is a solution of ()1.4j) . We 
compare it to another solution of the form 

(3.1) E = E + e a + e, n = n a + n 

where (e a ,n a ) is a solution of the homogeneous linearized equation (|2.1|) 
(|2.2p . with small boundary value on {z = 0} and such that n a and n are 
arbitrarily large at an arbitrarily small distance z > 0. Since we consider 



17 



periodic functions in x, it is sufficient to exhibit an example where x is one 
dimensional, which we now assume. We also assume, without restriction, 
that E is real. 

We choose (e a ,n a ) to be the solution of 



(3.2) 



id z e a + A x e a - En a = 0, 

{8 2 - A x )n a - A x {Ee a + E&) 



(3.3) 



n 



\t=0 



0, 



'|z=0 



a- 



sin(kt) 
J—' 



,ikx 



where the frequency k is large and a is a small parameter. Typically, we 
will choose a ~ k~ a so that the boundary data is of order k s ~ a in H s . Thus 



-ikx 



(3.4) e a {t, z, x) = u{t, z)e lkx + v(t, z)e 

(3.5) n a (t, z, x) = v(t, z)e ikx + ^zje~ ikx 

where (u a ,v a ,u a ) solve 



(3.6) 



id ; 



u 



k 2 u a - Ev a = 0, 



id z v a + k 2 v a + Ev a = 0, 
[ (d t 2 + k 2 )v a + k 2 E(u a + v a ) = 0, 



(3.7) ^ =0 = a^ =0 = 0, u* z=0 = a —^, vf z=o = 0. 

The solution (e a ,n a ) is amplified as described in Theorem 12.81 We con- 
struct correctors (e, n) such that (E, n) given by (|3.1|) is an exact solution 
of the nonlinear system, on a small domain Q where \h\ <C \n a \ and such 
that n a and thus n is arbitrarily large on a part of the boundary dVt. The 
equations for the remainders (e, n) in (|3.ip read 




i9 2 e + A x e - En = (e a + e)(n a + n), 

- A x )h - A x (Ee + Sg) = A x ((e a + e)(^ + 1)) . 



We add homogeneous initial-boundary data 

(3.9) h\ t=Q = d t h\ t=0 = 0, e\ z=0 = 0. 

so that E has the same boundary data as E + e a . 

We solve the nonlinear equation (|3.8p by Picard's iterations, in suitable 
Banach spaces of analytic functions which we now describe. 
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3.2 Construction of the correctors 

The approximate solution (e a ,n a ) has only two frequencies in x, +k and 
—k. The nonlinear interaction will create all the harmonics. This leads to 
consider functions of the form 



(3.10) u = J2 u p( t , 



z)e ikpx . 



Definition 3.1. Given parameters k > 2, 5 G]0, 1], b G]0, 1[ and s > 1, we 
denote by E the space of functions (|3.10p such that 

(3.11) |Up(t ' Z)| - C (lW 5<?>eWP 
with 



(3.12) (p) = max{2, \p\}, p = yj2k\E\ 2 zt, 
on the domain 

(3.13) O = {(t,z,x) G [0,1] X [0,6] x T : <5e p < l}. 
77ie 6esi constant in (I3.1ip is i/ie norm of u in E. 

Remark 3.2. The factor 5e p < 1 present for the frequencies ±/c in the ap- 
proximate solution is expected to intervene at the power \p\ in the harmonic 
pk by successive nonlinear interaction. Note that the series A3. lQf) define 
real analytic functions in x when (t, z) £ {5e p < 1}. 

The space E and the domain depend strongly on the parameters k, b 
and 5. However, to lighten notations, we do not mention this dependence 
explicitly. 

Lemma 3.3. Consider the solution (e a ,n a ) of (13. 2h ()3.3f) with a = jS. 

There is a constant Cq independent ofk >2,b G]0, 1], 7 > and 5 E]0, 1] 
such that 

i) f a = n a u a G E and h a = u a u a G E and 

(3-14) ||/ o ||e<C ^, H^He < c £ 

ii) for all v G E, u a v and n a v belong to E and 
(3.15) \\u a v\\ E < C j;\\v\\e, ||n a t!|| E < C 7||u||e- 
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Lemma 3.4. E is a Banach algebra and there is a constant Cq independent 
of k, b and 5 such that for all u and v in E: 

(3.16) \\uv\\ E < C ||u||e IMIe- 
Corollary 3.5. Introduce the notations 

(3.17) $(e,n) = (e a + e)(n a + n) 

(3.18) V(u,v) := (e a +u)(e^ + v). 

Then 

o 

(11 \ 

||$(e,n)|| E < C \^— + -\\n\\ E + 7||e|| E ) + ||e|| E ||n|| E 

2 

||*(«,t;)||E < °o(j^ + ?(II«I|k + IMIe) + IMIe ||i?||e 
Next, consider the linear problem 



(3.19) 



' id z e + A x e - En = f, 

{d 2 t - A x )n - A x (Ee + Be) = A x h, 
k n\ t=0 = d t n\ t=Q = 0, e\ z=Q = 0. 



Proposition 3.6. For f and h in E, the solution (e, n) of (13. 19ft belongs to 
E and 

(3.20) ||e|| E < dbWfU + dlnkWhU 

(3.21) ||n|| E < d |ipln(l M)||.f|| E + Cli Me 

where C\ is independent of k, b and 5 

Denote by (e,n) = T(f,h) the solution of (|3. 19|) . The equation ([3? 
reads 

(3.22) (e,h) =r(*(e,n),*(e,e)) :=.F(e,ra) := (£(e,n),J\f(e,h)). 
Proposition 3.7. For 7 > 1 and 

(3.23) 4CoC l7 2 (6+^) <1, 
(3-24) 4C 1 C 07 2 (^ + i)<|^| 2 , 
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T maps Bi x Bi into itself, where Mr denotes the ball of radius R in E, and 

k 

is a contraction for the norm 

1„ , 
IfIIe + tIp||e- 

i/7>2. 

Proof. By Corollary 13.51 for (e,n) G Bi x Bi and 7 > 1, there holds 

2 2 

||$(e,n)|| E <4Co^, ||¥(e, e)|| E < 4C p . 
Thus, by Proposition and (|333l) (133^) 

fc k 2 k 



\\£(e,n)\\ E <AC 1 C (b^ + ^)<l 



and 

..... ... 4CiCo7 2 /ln(M) x n 

||AA(e,n)|| E < —L^Lf-L^ + _) < 1. 



This shows that T maps B 1 x Bi into itself. 

Consider (e, n) and (e',n') inBi xBi. Denote by 5e = e— e', 5n = n — n', 
<5<3? = $(e,n) - $(e',n') etc. There holds 

||^|| E <2C 7(||<ye|| K + i||5n|| E ), 



|<5*|| E <^IHk 

k 



and thus, using Proposition 13.61 and (|3.23|) (|3.24|) . 

II^IIe < -(||<5e|| E + i||<5n|| E ), 
7 k 

l\\SAfh<^{\\Se\\ M + l\\Sn\W 

implying that T is a contraction if 7 > 2. □ 

Corollary 3.8. Assume that 7 > 2 and that (|3.23p and (|3.24p are satisfied. 
Then the equation (]3.8[) has a unique solution in Bi x Bj. 
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3.3 Choice of parameters, proof of Theorem 11.21 

Fix s = 1, e G]0, 1[, e' G]0, |e[ and cr > 0. We choose 

(3.25) 5 = k- a , b = k' e , ~f = k e '. 

Then the conditions (|3.23|) and (I3.24p are satisfied for k large and there is a 
solution 

E = E + e a + e, n = n a + n 
of the original system, with 

(3.26) l|e||E<p Ne<1. 
Therefore the first Fourier coefficient of n satisfies 

(3.27) MM) -n\{t,z)\ < 5 2 e 2p < 1 
while, by Theorem 12.81 for t > ^ and p large enough 

(3.28) K(M)| > np- 5 / 2 5e p . 

Note that 5e p = 1 for p = ^2k\E\ 2 zt = ahik thus for zt = ak^^nk) 2 
for some constant a > 0. For /c large, ak~ l (\n.k) 2 <C k~ e and therefore 

(3.29) T :=nn{5e p = 1} ^<D. 
On this set, 

he' 

(3.30) K(M)| >c'- - T . 

(In k) 2 

which tends to +oo with k. With H3.27I) . this finishes the proof of Theo- 
rem [L2j 

4 Proofs 
4.1 Lemmas 

We collect here several elementary estimates which will be used repeatedly 
in the sequel. 
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Lemma 4.1. There is a constant C such that for all A > 0, t > and all 

z > 0: 



z r t 



(4.1) / / e-^^dz'alt' < c\ ln(l + Azt), 

Jo Jo A 



(4.2) f f z'e-^dz'df 

Jo Jo A 



(4.3) / / t'e-^^dz'dt' < 

Jo Jo 



z ft 



Ct 



Proof. 



with 



Thus 



T e-^dt' = -^-H(Xz't) 
Jo A* 

H{u) = f e-^du' < C 
Jo 



l + u 



'* r^dz'dt! <C f — '— dz' = - r ''' <ls 



o jo jo 1 + Xz't X J 1 + s 

and (|4.1|) follows. The other two estimates are symmetric in t and z. Inte- 
grating in t' first as above we see that 



' f ' z'e-^dz'dt 1 <C I" tZ ' dz' < % [ Z dz' 
io Jo ~ Jo 1 + Xz't ~ X J 

implying (|4.2p . □ 
Lemma 4.2. For p > 1 and > 8, 

-ln(l + /ip 2 ) < ln(l + /A 
P 

Proof. The upper bound of the left hand side for real p > 1 is equal to c§^f\i 
if < 2/q with 7/o being the positive root of ln(l + t/q) = , where the 

function ln ^y V - reached its maximum equal to Co = • When \i > y 2 , 
the upper bound is attained at p = 1 and thus equal to ln(l + fj,). Because 
Vq < 8, the lemma follows. □ 



23 



4.2 Linear estimates 

In this section we prove the estimates of Proposition 13.61 Expanding the 
system (|3.19p in Fourier series in x and denoting by (e p ,n p , f p ) and h p the 
Fourier coefficients, yields the following system for U p := (u p ,v p ) = (e p ,eZp~) 



(4.4) 



' id z u p -k p u p - En p = f p , 
id z v p + k 2 p 2 v p + En p = g p , 
k (<9 t 2 + k 2 )n p + k 2 p 2 E(u p + v p ) = -k 2 p 2 h p , 



with g p = —f-p, plus initial boundary boundary conditions: 

(4.5) n p \ t=Q = 0, d t n p \ t=Q = 0, u p \ z=0 = v p \ z=0 = 0. 

Proposition 13.61 is an immediate corollary of the following estimates, where 
we use the notations 



(4.6) p(t,z) = yj2k\E\ 2 zt . 

Proposition 4.3. There is a constant C\ is independent of k > 2, b < 1 
and p such that for F p := (f p ,g p ) and h p satisfying on £1 

(4.7) \F p (t,z)\ < Ae {p)pit ' z \ \h p (t, z)\ < Be {p)p{t ' z \ 
the solution (U p ,n p ) of (14. 4|) ()4.5|) on £1 satisfies 

(4.8) \U p (t,z)\ <C x {bA + ^hx{hx8\)By^ z \ 

(4.9) \n p (t,z)\ < Ci(^ln(|ln<y|)A + ^(fc + lnfl ln S\))B)e^ z l 



Proof. We use Proposition 12.71 at the frequency kp. 

a) By (|2.32p . U p is the sum of three terms. The first one is bounded by 

(4.10) A [ Z e <PW>dz' < Aze {p)p(t ' z) . 

Jo 

When p = 0, only this term is present. When p ^ 0, the second term is 

U pfl : =£ *F P = [ [ S (t', z')F p {t -t',z- z')dz'dt'. 
JoJo 

Note the following identity for non negative real numbers: 
Vo7b> + VaJb" < y/(a' + a"){V + b"). 
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In particular, 

(4.11) p(t', z') + p(t -t',z- z') < p(t, z), 

and therefore, by Proposition 12.51 and f|4.Tj) : 



JoJo J(l + P') 



where ^ = v / 2kp\E\ 2 z't' and = 1 + *^ff z ' < 1 + /cp|£|V . 
When p > 2, 

z') - p' p = yj2(p 2 -p)k\E\ 2 z't' > yJkp 2 \E\ 2 z't' > ~pp{t', z'). 
When p = 1 , 

(l)p(i', z') - pi = 2p(i', z') - p(t', z') = p(t', z'). 
Thus, in any case, 

(4.12) {p)p{t\z')-p' p > l -pp(t\z'). 
Therefore 

\U p . (t,z)\ < Ae {p)p(t ' z) f K f\l + kplEfzy-^'^dz'dt'. 

JoJo 

Using Lemmas 14. II with A = ^kp 2 \E_\ 2 , yields 

(4.13) \U p , (t, z)\ < A(~ In (1 + Xzt) + kp\E\ 2 j)e™ < Aze^, 

where we have used that j In (l + Xzt) < zt. 
Similarly U Pi i = E\ * h p satisfies 

\U Ptl (t,z)\<Be^ rfe-PP'JM^^dz'dt' 
JoJo J(l + p' p ) 



Zf-t 



<Be^ p / / kpe-2 pp 'dz'dt' 
JoJo 

<B^ln(l + iA|£| 2 .t) e®'. 



On \k\E\ 2 zt < p 2 < | ln5\ 2 . Thus, for 5 > e 3 , Lemma S3] implies that 
(4.14) \U p!l (t,z)\ < ln(l + |ln5|V P>P < j^MIM) e®'. 
With KlOh and (H~T3l) . this implies the estimate (fl~8j) . 
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b ) Similarly, when p = the estimate for no is immediate. When 
p ^ 0, by (|2.33p . n p is the sum of two terms. Up to constant factors, the 
first one is a convolution of F p by £\ (computed at the frequency kp). Using 
again (|4.11j) and (|4.12p this term satisfies 



\n Pil (t,z)\<Ae^ / / Kpe . dz'dt' 



(4.15) 



070 yj^+pfp) 



<Ae {p)p [ [ kpe'^' dz'dt' < A-^- ln(| In 5\)e^ p . 



JOJO 

The second term in n p is the convolution of h p with £2 and satisfies 

\n p>2 (t,z)\<Be^) / / e-^ k p dz'dt' 

J0J0 f' pS J{l + p' p ) 

with jr < (f + 5^). Thus, by Lemma gj] 



(4.16) 



\np >2 (t,z)\ < Be^ p [ ( k 2 p 2 (^-+t')e-^ t '' z 'Uz'dt' 



kp 

kt 1 

^'(r^ + jlj5 MIM >)'"'' 



□ 



5 Analytic solutions of the Goursat problem 



Consider 




id z E + A X E = nE, 
(d? - A x )n = A X \E 



on t £ [0,T], z £ [0, Z\ and x G T 2 , together with initial boundary conditions 

(5.2) ri| i=0 = 0, d t n\ t=0 = 0, E z=0 = E . 

We prove that this Goursat problem is locally well posed in spaces of 
analytic functions, following the general approach presented in |llj . One of 
our objective is to give an explicit lower bound for the domain of existence, 
depending on the norms of the boundary data (see Remark I5.6I below) . 



26 



5.1 Spaces 

Fix s > and equip H S (T) with a norm || • || s such that 

(5.3) lis < \\ u \\s \\v\\s- 

Definition 5.1. Given a (formal) power series 4>(x) = ^ </> n x n in one vari- 
able x, with nonnegative coefficients, we say that u(x) <C <f> if for all a £ Z 2 , 
[|9£it|| a < |«| ! 0|a|- 

Lemma 5.2. i) If u <C 4>, then d Xj u <C <ft' ■ 
i) If u <C 4> and v <C ijj, then uv <C 4>ijj, 

We will apply this definition to functions u and power series which also 
depend on t and z, seen as parameters. In particular, given a power series 
</>(x), we will consider power series 

(5.4) 0a(*,x) = <Xx + Az), 

and use the notations z, x) <C 0a(<z> x ) if u(i, z, •) <C <f>\(z, •). 
Lemma 5.3. Ifu(t,z,x) <C 4^r\{ z i x ) 

u(t,z,x)dz « -cf) R ,\(t,z,x). 

In particular, we consider the power series (see |llj ) 
~ x n 

(5.5) y?(x) = co V 2 - , ^h,a(^x) = <p(Rx + Xz). 

n=0 

where Co is such that y? 2 <C We also introduce the notations 

oo 

(5.6) <p(Rx + Xz) = (fin(Xz)R n x n . 

n=0 

Of course the explicit value of <p n can be deduced from (|5.5p . We also note 
that 

oo 

(5.7) *4,a(z,x) = J> + 1 W(Az)iT+V. 

n=0 

and by integration of (p'(Rx + Xz) in z, we see that 

r z i i 

(5.8) / ip n+1 {Xz)dz = - — — — (ip n (Xz) - <p n (0)) < - — — — ip n {Xz) 
Jo (n + l)A v (n + l)A 
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Definition 5.4. Denote by Er^ the space of functions u on [0, T] x [0, Z] x 

T 2 , with Z = j, such that 

(5.9) u(t,z,x) <C C(p R) \(z,x) 

/or some constant C > 0. 

Similarly, ¥r\ denotes the space of functions u on [0,T] x [0, if] x T 2 , 
sitc/i that 

(5.10) u(t,z,x) <C C^ )A (2,x) 
/or some constant C > 0. 

The norms in E and F are the best constant C in (15. 9p and (15.101) . 
For the boundary data Eo(t, x) on [0, T] x T 2 , we consider the spaces E# 
of functions u(t,x) <C C(p(Rx). 

Theorem 5.5. Lei i?o £ E_r. Then, the problem (|5.1j) /ias a unique solution 
(E,n) G E/j a x F/j a if X> CRT\\Eq\\^ where C is a constant independent 
of T , R and eo . 

Remark 5.6. Functions in E^ are defined for i2|Imx| + X\z\ < 1. Thus 
measures the width of the complex domain where the boundary data is 
defined and A -1 is the order of the length of propagation in z. In particular, 
for boundary data Eo(t,x) = e(t)e tkx then, for R « k, \\Eq\\e r < E = 
||e||ioo. Theorem 15.51 asserts that the length of stability Z = A -1 satisfies 

(5.11) ZTkE 2 < 1, 

which is very similar to the condition p < 1 that was used in Section 3. 

5.2 Resolution of (EU) in E 

Consider the wave equation 

(5.12) d\n - A x n = A x h, n\ t=Q = d t n\ t=0 = 0. 

Lemma 5.7. For h G ^r,\ , the solution n of (|5.12p belongs to F^a and 
(5-13) ||n|| F ^ < C T\\h\\ ERX 

where Cq is independent of T, R and A. 
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Proof. The mean value of n (the 0-th Fourier coefficient) vanishes and there- 
fore the H s norm of n is equivalent to the i/ s_1 norm of d x n. Thus, by 
standard energy estimates for the wave equation, there holds 

\\d«n(t,z,-)\\ s <C f \\d«h{t\z,-)\\ s+1 dt'. 
Jo 

Thus, for \a\ = n, 

\\d%n(t,z,-)\\ s <CT\\h\k RtX (n+1)! R n+1 cp n+1 (Xz). 

Therefore 

n(t,z,x) < Cr||/t||]E K:A ip' RX (z,x) 
and the lemma follows. □ 

We consider next the Schrodinger equation 

(5.14) id t e - A x e = f, e\ z=0 = e . 

Lemma 5.8. For f G ^r,\, and eo G E_r, the solution e of (|5.14p belongs 
to Mrx and 

R 

(5-15) IMIeh.a ^ yll/lkii,A + ll e ol|E_R- 

Proof. Standard energy estimates imply that 

||^e(t s 0, Oil, < f 011^' + ll^eo(t, -)ll. 



Thus, for \a\ = n, 

\\d^e(t,z,-)\\ s <\\h\\ ¥RA (n + l)lR n+1 I <p n+1 (\z')dz' 



R n+l 
Therefore 

e{t,z,x) < yH/HF^ ipR,\(z,x) 
and the lemma follows. 

□ 

Lemma 5.9. For e\ G Er^, 62 G E_r 5 a an <^ n 6 F^, there holds e±e2 G 
Fr,a and ne\ G F^a 

||eie 2 ||E Je , A < ||ei||E fl , A ||e2||E_ R , A , 
||nei|| F p A < ||ei||E fliA INK,*- 
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5.3 Proof of Theorem 1531 

Denote by M(h) the solution of (|5.12|) and by So(eo) + <S(f) the solution of 
(|5.14p . To prove Theorem 15.51 it is sufficient to show that there is C > 0, 
such that for A > CiJTHE'oll^, the equation 

(5.16) E = S(E ) + T (E, E, E) 

has a unique solution in IEr^, where T is the trilinear operator 

T{u, v,w) = S (uAf (vw)) . 

The estimates above show that T maps (E^a) 3 to Ej^a and that 

Co RT 

\\T(u,v,w)\\e rx < -^—\\u\\e Ri Jv\\e r< Jw\\e r<x . 
From here, standard Picard's iterates imply Theorem 15.51 

6 Remarks and comments 
6.1 Nonconstant backgrounds 

The boundary data we have considered in Sections 2 and 3 are of the form 
(6.1) l{t>o}(# + e(t,z)) 5 

thus have a jump at t = 0, which is not physical. Note that arbitrary 
functions E(t) together with n = are solutions of (jl.4p . One could take 
them as background state. We now briefly sketch how one can start the 
analysis. Again, there is no restriction in assuming that E is real. The 
linearized systems (|2.3p and (|2.4p are unchanged, except that E_ is now a 
function of time. Therefore, we now proceed by using the Fourier-Laplace 
in z rather than in time: we_extend u, v , n by for z < 0. Note that the 
extension u satisfies d z u = d z u + u$5 z= q. Therefore, the Fourier-Laplace 
transforms satisfy 



(6.2) 



— (( + k )u — En = iuo, 
(-( + k 2 )v + En = ivo, 
{ {d? + k 2 )h + k 2 E(u + v) 
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Hence, 

(6.3) 

and 

(6.4) 

with 

(6.5) 
(6.6) 



f En iuQ 

u = — 



C + k 2 ( + k 2 ' 
En ivo 



(-k 2 Q-k 2 ' 



(df + r(l + ap{t))n = k 2 h, 



a = \E\ 



2k 2 



C 2 -k^ 



Pit) = \E(t)\ : 



Denote by M(t, s, a) the fundamental solution of the second order o.d.e. 
(6.7) (d 2 + k 2 (l + ap))N = 0, A/] t=s = 0, d t M\ t=s = l. 
Thus 



(6.8) 



n(t) = k 2 I U(t,s,a)h(s)ds. 
Jo 



Performing the inverse Laplace transform, we get 

(6.9) n= £\(t, s, z)uo(s)ds + / £2(^7 s, z)vo(s)ds 

Jo Jo 

with 

(6.10) £i(t,s,z) 



2vr 



(6.11) 



£ 2 (t,s,z) 



iEjs) 
2vr 



<U(t, 



s,a, 



k 2 d( 
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Lemma 6.1. N is an entire function of a. Moreover, for \a\ small and 
< s < t < 1 there holds 

(6.12) Af(t, s,a) = l (a(t, s, a)e ikv &'*) _ p(t, s, a) e -^ (M ' a) ) 

with 



(6.13) d t (p = y/TTap, <p\ t=s = 0, 

and a and (3 are 0(1) as well as their time derivatives. 

Because \a\ — > Oas |£| — > oo, U is holomorphic in (" and bounded for 
large £. Thus, one can close the integration path in (|6.10jl and (|6.1ip : for 
instance, for z > 0: 

1 f izCrrr, _ ^ ^< 



(6.14) Si(t,s,z) = ^- [ e izC U(t,s,a 

2ir J r 



(-k 2 



where T is a large circle in the complex plane. 

From here, using the Lemma above, one can repeat most of the com- 
putations performed in Section 2, with only technical additional difficulties. 
The important point is that the amplification is now driven by the kernel e p 
with 

(6.15) p:= ^J2k{z - z')(P(t) - P(t% P(t)= [ \E(s)\ 2 ds. 

Jo 

Of course, when E is constant we recover the previous results. On an interval 
where \E_\ > 0, there holds P(t) — P(t') w (t — t') and all the estimates of 
Section 2 are unchanged. 

In this direction, another interpretation of the constant case E, different 
from (|6.ip is that the boundary data is 



(6.16) E(t) + e(t,x)l {tm 

where E vanishes for t < — Tq < and is equal to the constant E for 
t > 0. This means that the perturbation e starts when the background is 
stabilized at a constant value. This is not physically realistic either, but 
the expectation is that this analysis is qualitatively relevant for physical 
interpretations. 
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6.2 The amplification rate 

We come back to the case where E is constant. The solution of the ho- 
mogeneous linearized equations (|2,4p is given in Lemma 12.11 It involves 
convolution by functions defined by contour integrals (|2.15|) . These inte- 
grals were estimated in Section 2 using deformations of contours and the 
saddle point method. According to Lemma 12.41 a model for the integral 8i, 
is 

(6.17) = J- *C X = k+ '- |2 



2iir J r C-fc 2(C - k) 

Note that 

(6.18) T\ = e l(tk - zk2) J(p), p = y / 2k\E\ 2 zt 

where J is the bessel function 

n=0 



2*xJr C Jo ^nl 2 



This is indeed the core of the analysis performed in Section 2. 

Following this path, a model for the solutions of the homogeneous equa- 
tion (El) is 



(6.19) u(t,z)= e l{{t - tr)k ~ zk ' 2) J{^2k\E\ 2 z{t - t'))u {t')dt' . 



The kernel has an amplification part given by J but it also has an oscillating 
part given by e^ i_i ^ k . In Theorem 12.81 we have chosen uq such that it has 
itself an oscillation e %kt to eliminate the oscillations in (|6.19p and maximize 
the amplification. For general uq the actual amplification results from a 
delicate balance between J and the oscillations. Typically, one integration 
by part in the integral above, wins a factor k~z and k~ l times a derivative 
of no- Thus the amplification e p is cut down by a factor k~ m where m 
is the smoothness of Uq. This can also be seen on the Fourier Laplace 
representation of u: 



e i(tk-zk 2 ) r »»LglV d( 



(6.20) u(t, z) = — jf e i[ *-— h (k + C) 

In any case, there is one clear conclusion : the more amplified uq are those 
who have a nontrivial oscillation at frequency k. 
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This can also be seen in a different way. Consider the Fourier-Laplace 
expression of the solution: 



(6.21) U(z,() = e izA m (O- 

For £ = (r — 27) ^ {—k, +k}, the eigenvalues of A are 



(6.22) M (C) = ±fcVl + 



6' ISP ' 



Consider real time frequencies C, = t ^ {— k, +k}. When 1 + i > 0, the 
eigenvalues are real and there is no amplification. When 1 + | < 0, that is 
when 



(6.23) < k 2 - t 2 < \E 



2 



the eigenvalue is purely imaginary and this suggests that the amplification 
factor for oscillations with frequency r is 



(6.24) e 2 l Im ^l, |Im M | = ^.^^ 



T 



2 



In particular, only real frequencies close to k are amplified. 

Of course, considering only real frequencies is not consistent with the 
forward evolution problem under consideration and the condition that Uo 
vanishes in the past. e lzA ^ has essential singularities at +k and —k, and 
one has to consider integration on complex path as in Section 2, to turn 
around the singular points. Then the imaginary part of the phase to consider 
is tlm £ — zlm \x. In Section 2 we have chosen optimal contours in order 
to obtain sharp estimates. In particular, the amplification (16.190 is not 
correct. But the computation suggests that the amplification is mainly due 
to frequencies close to k, and this is correct. 



6.3 The physical context and physical values 

The Zakharov's equations [12] have been introduced at the beginning of the 
70's, to describe electronic plasma waves. They couple the slowly varying 
envelope of the electric field and the low-frequency variation of the density 
of the ions. When modeling the propagation of a laser beam in a plasma, 
several phenomena occur. One has to take into account the laser beam 
itself, the Raman component and the electronic plasma waves (see [2j [3] for 
example). For laser propagation or for the Raman component, one often 
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uses the paraxial approximation and the Zakharov system that couples the 
envelop of the vector potential A of the electromagnetic field A to the low- 
frequency variation of the density Sn of the ions reads 



(6.25) 



LO, 



pe 



L0 

(df-c 2 s A x )Sn 



2to 



2n uj 



■5nA, 



LO 



pe 



47rm,c 2 



A x \Af 



its wave number and to pe the 



where loq is the frequency of the laser, ko 

plasma electronic frequency; they are linked by the dispersion relation to 2 = 
u r>e + ^o c2 where c is the speed of light; no the mean density of the plasma, 



J pe 

rrii is the mass of the ions and c s the sound velocity in the plasma, 
suitable units for no, the plasma frequency to pe is given by 



In 



(6.26) 



UJ- 



pe 



4-7re 2 no 



where — e is the charge of the electron and m e its mass. 

As above, z is the space component in the direction of propagation of 
the laser beam and x denotes the space components in directions that are 
transversal to the propagation. 

Introducing the dimensionless quantities 

A 



(6.27) 



n 



5n/n Q , 



a = —r 



A- 4 



nipcc. 



and using (|6.26p . the system reads 



(6.28) 



,w 



1 



1 



d t + k d z )a + -A x a 



to 



pe 



2c 2 



n a, 



{-^df - A x )n 



m f 



AJa\ 



rrn 



In addition to the pulsation ojq and the wave length Ao = character- 
istic quantities for the laser beam are its duration To, its transversal width 
i?o (supposed to be smaller than the transversal dimension of the plasma), 
the expected length of propagation Zq and also the characteristic dimension 
of the transversal variations (speckles) A which correspond to variations of 
the spectrum by frequencies k ~ 2j£ In the application we have in mind 
A -C i?o and both for theoretical and computational reasons (spectral meth- 
ods) it makes sense to assume periodicity in x with period X satisfying 
A < X < #0- 
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With this data in mind, rescale the space-time variables introducing the 
characteristic transversal width X: 



x 



t= x 



{df - A £ )n = ^A £ \a\ 2 , 



2uj c s X 



LV pe X 



(6 - 29 > — *• -2W" 

In these variables, the system reads 

i{edi + dz)a + A^a = a 2 na, 

(6.30) I _ . x m P 

with 

(6.31) e = — — , a 
With 

(6.32) n = a 2 n, a = a 

we obtain the dimensionless system: 

f z(e<9 f + <9 5 )a + A s a = no, 
K '" ,) \(3$-A £ )h = A,\d\ 2 . 

For this system, the typical amplification factor is 




(6.34) 



2\a\ 2 kfo Zq 



where k is the scaled frequency of speckles, tq the scaled time of propagation 
under consideration and Zq the scaled length of propagation. Scaling back 
to the original variables, there holds 



(6.35) 



thus 



2k X 2 ' 



k = Xk, f 



c s Tq 
X ' 



Proposition 6.2. The dimensioned amplification factor p of Section 2 is 
given by 



(6.36) 



2 _ \A\ 2 m e c s uj 

P — 79 < 

A nijC 



T tqZq. 
k 
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For the physical significance of the analysis it is important to evaluate the 
various constants. The velocity of light is c = 3 10 8 ms _1 , the sound velocity 
of the electrons is of order c s « 0.005c = 1.5 10 6 ms~ 1 . With e = 1.6 10 _19 C 
and m e = 0.9 10 _30 Kg, the scaling factor for A is 

(6.37) A=^£«2.510 3 V. 



Typical values of io pe are of order 10 15 s 1 and 2k « 10 4 . Thus 

; 2 
rrijC 2 



(6 .3 8) 2 ^W^ ^ 4 iq15 m _ lg _ 1 



A typical value of Ao is Ao = 0.35/iin corresponding to a wave number 
ko = ^ = 1.8 10 7 m _1 . The pulsation is given by the dispersion relation 
lOq = ujp e + A;qC 2 , yielding w = 5.5 10 15 s~ 1 . 

By construction, the relative variation n = 5n/riQ is small compared 
to 1 and the consistency of (|6.28p requires that 7^|a| 2 must also be small 
compared to 1, that is \a\ <C 100 yielding the upper bound for A 

(6.39) \A\ < 10 5 V. 

The envelop of the electromagnetic field is linked to A by the polarization 
relation 

(6.40) E = i—A 

c 

yielding the bound 

(6.41) \E\ < 10 12 Vm" 1 , 

which allows for very intense fields. 

Typical values for the dimension of the speckles is A ~ 10/um, corre- 
sponding to frequencies k = 2j£ ps 6 10 5 m _1 . A constraint for the validity of 
the paraxial approximation which sustains the Schrodinger envelop equation 
in (f6T25j) is that 

(6.42) kZ < k R 

With ko/k ~ 30 as above and a diameter Rq ranging from 1mm to 1cm, this 
allows for propagations along distances Zq of order 1 to several cm. In this 
range of k, we finally obtain: 

(6.43) p= \^-y/pT Zo with /3«10 14 m^^" 1 
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The duration of the laser impulsion tq may be of order 10 -8 to 10~ 10 s. 
Note that for a 1cm long propagation of the laser beam, the time elapsed is 
of order 310 _11 s. For Zq of order 1cm and tq ranging from 10~ 8 to 10~ 10 s, 
we obtain that 

(6.44) p = 7 ^ 

with 7 ranging from 10 to 100. 

Thus, for p to be small, one must have 

(6.45) |A|<10V, or \E\ < 10 8 Vm" 1 , 

which does nor correspond to high intensity beams. 
Note also that with X ~ 50 pm, there holds 

(6.46) e ss 9, f « 3. 

In conclusion, we see that the values obtained in f|6.45[> are compatible 
with the paraxial approximation and the derivation of the model (|6.25|) . 
However, they do not allow for high intensity or high energy beams. This 
suggests that the model could be unadapted to such situations. 
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